Skyrmions with vector mesons in the hidden local symmetry approach 
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The roles of light p and uj vector mesons in the Skyrmion are investigated in a chiral Lagrangian 
derived from the hidden local symmetry (HLS) up to 0(p 4 ) including the homogeneous Wess-Zumino 
(hWZ) terms. We write a general "master formula" that allows us to determine the parameters of the 
HLS Lagrangian from a class of holographic QCD models valid at large N c and A ('t Hooft constant) 
limit by integrating out the infinite towers of vector and axial-vector mesons other than the lowest 
p and uj mesons. Within this approach we find that the physical properties of the Skyrmion as 
the solitonic description of baryons are independent of the HLS parameter a. Therefore the only 
parameters of the model are the pion decay constant and the vector meson mass. Once determined 
in the meson sector, we have a totally parameter-free theory that allows us to study unequivocally 
the role of light vector mesons in the Skyrmion structure. We find, as suggested by Sutcliffe, that 
inclusion of the p meson reduces the soliton mass, which makes the Skyrmion come closer to the 
Bogomol'nyi-Prasad-Sommerfield (BPS) soliton, but the role of the uj meson is found to increase 
the soliton mass. In a stark contrast, the A-N mass difference, which is determined by the moment 
of inertia in adiabatic collective quantization of the Skyrmion, is increased by the p vector meson, 
while it is reduced by the inclusion of the uj meson. All these observations show the importance of 
the uj meson in the properties of the nucleon and nuclear matter in the Skyrme model. 

PACS numbers: 12.39.Dc, 12.39.Fe, 14.20.-c 



I. INTRODUCTION 

In accessing dense baryonic matter, one possible ap- 
proach that unifies both the elementary baryons and 
multi-baryons system was proposed in Refs. [HQ ■ hi this 
approach, starting with a chiral Lagrangian, the single 
baryon is generated as a Skyrmion and multi-Skyrmions 
are put on crystal lattice to simulate many-baryon sys- 
tems and dense matter. 

However, the previous works in this approach suffer 
from the modeling of the effective chiral Lagrangian and 
determination of the low energy constants (LECs). This 
problem becomes serious when one considers higher or- 
der chiral Lagrangian or introduce more mesonic de- 
grees of freedom, which prevents systematic studies for 
Skyrmion properties. On the other hand, when one starts 
with holographic QCD (hQCD) models and integrates 
out infinite towers of mesons such as vector and axial- 
vector mesons except a few low-lying mesons as done in 
Refs. [1,01: h leads to a chiral Lagrangian with the val- 
ues of all the LECs fixed by only a few phenomenological 
inputs. Therefore, the number of parameters is dras- 
tically reduced, which allows the studies on Skyrmions 
in a systematic way. In the present article, we use the 
general master formula derived in this way to determine 
the LECs of the 0(p A ) terms of the hidden local sym- 
metry (HLS) Lagrangian, which leaves undetermined by 
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theory the pion decay constant, vector meson mass, and 
the HLS parameter a. As will be shown, all physical ob- 
servables that we will consider are independent of a. So 
by fixing the pion decay constant and the vector meson 
mass from the meson sector, we are able to study the 
Skyrmion properties in a totally parameter-independent 
way, a feat that as far as we know, has not been achieved 
before. In addition, we work with the chiral Lagrangian 
including the homogeneous Wess-Zumino (hWZ) terms 
for studying the role of the uj vector meson in the prop- 
erties of a single Skyrmion. This is our first step towards 
more complete studies on Skyrmions for nucleon struc- 
ture and baryonic matter. The main results of this work 
were quoted in Ref. Q, and here we provide the details 
of the calculations in this model as well as more detailed 
analyses on the Skyrmion mass and size. 

The Skyrme model @, 0] is the nonlinear sigma model 
stabilized by the Skyrme term, a four-derivative term, 
where the baryons emerge as stable field configurations 
with a non-trivial geometrical structure. It is well ac- 
cepted that the nonlinear sigma model captures the 
physics of QCD at very low energy scales and as the en- 
ergy scales increase vector mesons should be excited. An 
elegant way to describe the vector meson physics is the 
Hidden Local Symmetry (HLS) in which the vector 

mesons emerge as the gauge bosons of the HLS. Further- 
more, as the energy scale goes up, infinite number of local 
symmetries appear and the corresponding gauge fields 
are identified with the infinite vector and axial-vector 
mesons. These infinite number of hidden gauge vector 
fields together with the pion field in 4-dimension (4D) can 
be dimcnsionally de-constructed to 5-dimensional (5D) 
Yang-Mills (YM) action in curved space with the 
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5th dimension being the energy scale. 

A similar situation exits in the gravity sector (that 
is referred to as "bulk" sector) of gravity /gauge (holo- 
graphic) duality that comes from string theory. (See, 
for example, a recent review addressed to hadron physi- 
cists in Ref. [l2[ and references therein.) When Kaluza- 
Klein(KK)-decomposed to 4D and all the KK modes ex- 
cept the lowest lying vector mesons and pscudoscalar 
mesons are integrated out in a way consistent with hidden 
local symmetry from the bulk sector, this model is dual 
to the HLS up to 0(p 4 ) @, 13 . This dual (bulk-sector) 
model is justified in the large number of colors N c , large 
't Hooft coupling A = Qym-^c limit and the chiral limit 
where the quark masses vanish. Furthermore, baryons 
can be described as solitons in the holographic QCD (]~3l — 
Il7j |. The hQCD model includes only two parameters, A 
and the KK mass Mkk, which can be fixed from meson 
physics. This, therefore, enables a parameter-free cal- 
culation of the Skyrmion properties with vector mesons 
and provides a way to perform a systematic study on the 
role of vector mesons in the Skyrmion structure. Here, 
the five-dimensional Chern-Simons (CS) term that is re- 
sponsible for the anomalous part of the HLS Lagrangian 
is a topological quantity and, therefore, is free from the 
warping of the space-time. As we shall see below, the CS 
term is very important to understand the role of the to 
meson in the soliton structure. 

In the literature, Skyrmion has been studied based on 
the 0(p 2 ) HLS Lagrangian as in Refs fll^. This model 
has three parameter, f w , g, and a, where /„■ is the pion 
decay constant, g is the HLS gauge coupling constant, 
and a is a free parameter in the HLS. (See Section [TT1 for 
the definition of these parameters.) The HLS parame- 
ter a is normally taken to be 1 < a < 2 [§4lol|. In free 
space a ~ 2 is preferred but in hadronic medium at high 
temperature and/or density, one gets a ~ 1 |l£j. The 
dependence of a on circumstances hinders systematic in- 
vestigation on the properties of a single Skyrmion and 
baryonic matter. For example, the soliton mass reported 
in Ref. (l9j within a p-meson stabilized model is 

M so i = (667 - 1575) MeV (1) 

for 1 < a < 4 with m T = 0, and the pion mass effect is 
found to be small in the soliton mass. This shows that the 
ambiguity in the value of a results in a large uncertainty 
in the soliton mass. 

Furthermore, the description of baryons as Skyrmions 
is supported by the large N c limit [221 . I n the HLS, the 
higher order terms such as the 0(p 4 ) terms are at 0(N C ) 
like the 0{p 2 ) terms0 As a result, in the N c counting, 
these higher order terms should be taken into account. 
However, including the higher order terms inevitably 
calls for more complicated form of the Lagrangian and 



The loop corrections from the 0(p 2 ) Lagrangian to the 0(p 4 ) 
terms arc 0(N®) and therefore are sub-dominant. 



uncontrollably large number of low energy constants. In 
this paper, these constants will be determined in a con- 
trollable way by using a master formula. 

This paper is organized as follows. In Sec. |Hl we in- 
troduce the HLS Lagrangian up to 0(p 4 ) including all 
the hWZ terms. The soliton wave functions are con- 
structed and the collective quantization method is also 
briefly explained. We show a general master formula to 
determine the parameters of the HLS Lagrangian induced 
from a class of 5D gauge models including hQCD mod- 
els in Sec. Mil In Sec. IIV| we present our results on the 
Skyrmion mass and size as well as the moment of in- 
ertia calculated in the present work. Here, we consider 
two models for the parameters, which include the HLS 
induced from the Sakai-Sugimoto (SS) model and the 
HLS induced from the Bogomornyi-Prasad-Sommcrficld 
(BPS) model. The results from these two parameter sets 
are discussed and compared. Section |V] contains a sum- 
mary and discussion. The complete explicit expression 
for the soliton mass and the equations of motion for the 
static fields are given in Appendix [A] The moment of 
inertia and the associated equations of motion for the 
excited fields are collected in Appendix [Bl 



II. SKYRMIONS FROM THE HIDDEN LOCAL 
SYMMETRY 

In order to study the role of vector mesons in 
Skyrmions, we first briefly introduce the chiral effective 
Lagrangian with vector mesons referring for the details 
to Refs. [l(| H3l- Here, we consider both the iso-scalar 
oj meson and the iso-vector p meson as well as the chiral 
field as explicit degrees of freedom in the theory. These 
vector mesons are introduced as th e gaug e bosons of the 
HLS of the nonlinear sigma model j8l4l0j . 

The full symmetry group of our effective Lagrangian 
is Gf u ii = [SU(2) L x SU(2) R } chiral x [C/(2)]hls with 
[E/(2)]hls being the HLS. In the absence of the exter- 
nal sources, the HLS Lagrangian can be constructed by 
making use of the two 1-forms a\\^ and a± fl defined by 

«±p = ki ( D "Z*& - > ( 2 ) 

with the chiral fields £ L and which are written in the 
unitary gauge as 

CI = ^ : = € = e" /2/ " with it = tvt, (3) 

where t are the Pauli matrices. The covariant derivative 
is defined as 

D^r.l = (d„ - iV^ R ,L (4) 

with Vfj, being the gauge boson of the HLS. This is the 
way to introduce vector mesons in the HLS, where the 
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vector meson field is 



with 



Pn = Pti 



P° 



(5) 



(6) 



Then one can construct the chiral Lagrangian up to 
0(p 4 ) as 



AiLS = £(2) + £(4) + £anom, 

which is our working Lagrangian. Here, 



(7) 



£ (2) - tfTr (a ± ^ ± ) + af^Tr [a^ j - ^Tr (^V'"') 



1 



(8) 



where / ff is the pion decay constant, a is the parameter of the HLS, g is the vector meson coupling constant, and the 
field-strength tensor of the vector meson is 



V, 



d^-d^-iiV^Vv}. 



(9) 



In the most general form of the 0(p i ) Lagrangian there are several terms that include two traces in the flavor space 
such as the 2/io~2/i8 terms listed in Ref. [lO]0 These terms are suppressed by iV c compared to the other terms in the 
Lagrangian and are not considered in the present work. Then the 0(p 4 ) Lagrangian which we study in this paper is 
given by 



£( 4 ) = C 



(4)y 



C 



(4), 



(10) 



where 



C 



(4)y 



C 



(4)z 



2/i Tr 
+ %Tr 
+ y 8 {Tr 

V~LL * 1/ 



2/2 Tr 



2/6 Tr 

+ Tr 



*_L At u: _L I ,<- 



iz 4 Tr 



iz 5 Tr 



V 1 A u 



+ y 7 Tr 
} + % Tr 



y 4 Tr 



(11) 

(12) 



r 



In the present work, we also consider the anomalous par- 
ity hWZ terms that are written as 



N c 

167T 



where 



Ci = iTr[&l& R - d R d L ], 

C 2 = iTr[d L d R d L d R ], 

C 3 = Tt[F v (d L d R - d R d L )] . 

in the 1-form notation with 



(13) 



(14a) 
(14b) 
(14c) 



In order to study the properties of the soliton obtained 
from the Lagrangian (J7J, we take the standard parame- 
terization for the soliton configuration. For the pion field, 
we use the standard hedgehog configuration, 



£(r) = cxp 



F(r) 



(16) 



The configuration of the vector mesons are written as [2C 

G(r) 



gr 



(rxr). (17) 



olt. = an — a, 



tX R = OK || -|- Uj_, 

F v =dV - iV 2 



For the baryon number B = 1 solution, these wave func- 
tions satisfy the following boundary conditions: 



(15) 



F(0) = Tr, F(oo) = 0, 
G(0) = -2, G(oo) = 0, 
W'(0) = 0, W(oo) = 0. 



(18) 



2 Another example of this kind is Tr [Ay 1 ] Tr [ay J that generates 
the mass difference between the p and ui mesons. 



Given the Lagrangian and the wave functions, it is now 
straightforward to derive the soliton mass M so \. The 
explicit expression for the soliton mass is given in Ap- 
pendix |XJ Minimizing the soliton mass then gives the 
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coupled equations of motion for the wave functions F(r), 
W(r), and G(r). These are also given in Appendix 1X1 

The classical configuration of the soliton obtained 
above should be quantized to describe physical baryons 
of definite spin and isospin. Here, we follow the standard 
collective quantization method [24j , which transforms the 
chiral field and the vector meson field as 



V fl (r)^V^r,t)=A(t)V^r)A^t), 



(19) 



where A(t) is a time-dependent SU(2) matrix. We de- 
fine the angular velocity 1~2 of the collective coordinate 
rotation as 



ir-n = A\t)d A(t). 



(20) 



Under the rotation (TT9"]) , the space component of the w 
field and the time component of the p field, i.e., uf and 
p°, get excited. The most general forms for the vector- 
meson excitations are written as |20j 



P a (r , t) = A{t) - [t • « & (r) + <? • r n • r & (r)] A* (t), 
9 



w *(r,t) = ^(nxr) i , 



(21) 



With these wave functions the moment of inertia can 
be calculated and its explicit expression is given in Ap- 
pendix [Bj It is then straightforward to obtain the Euler- 
Lagrange equations for the wave functions, £i( r )) ^2( r )i 
and v?( r ) by minimizing the moment of inertia, and the 
results are also given in Appendix [SJ The boundary con- 
ditions imposed on the excited fields arc 



£(o) 

V(0) 



£2(0°) 



= 0, 

0, 



(22) 



and £i(r) and £ 2 ( r ) at r = satisfy the constraint, 

2£i(0)+C 2 (0)=2. (23) 

In the adiabatic collective quantization scheme, the 
baryon mass is given by 

M = M 8oI + ^ = M sol + M±}± (24) 



21 



21 



where i and j are isospin and spin of the baryon. Then 
the A-A^ mass difference reads 



Hi 



M A -M N = —. 



(25) 



The baryonic size of a baryon should be computed by 
the baryon number current of the Skyrmion. However, in 
order to intuitively see the effects of the vector mesons on 
the Skyrmion size in a simple way, here we consider the 
winding number and energy root mean square radii. The 



root mean square (rms) radius of the winding number 
current is defined by 



(r 



2\l/2 
IW 



d 6 rr z B u (r) 



1/2 



(26) 



where B° (r) is the time component of the winding num- 
ber current that is explicitly written as 



B° 



1 



27t 2 r 2 



F' sin 2 F. 



(27) 



We define the energy root mean square radius (r 2 ) 1 ^ 2 as 



(r 



2\l/2 



M so l 



-1 1/2 



d s rr 2 M soi {r) 



(28) 



where M so \(r) is the soliton mass (energy) density given 
in Appendix [A"l 



III. HIDDEN LOCAL SYMMETRY INDUCED 
FROM HOLOGRAPHIC QCD 

A. Master formula 

In this section, following Refs. [1,01 , we provide a gen- 
eral master formula to determine the parameters of the 
HLS Lagrangian by integrating out the infinite towers of 
vector and axial- vector mesons in a class of hQCD models 
expressed by the following general 5D action: 



oDBI 



s; 



cs 



(29) 



where the 5D Dirac-Born-Infeld (DBI) part Sf BI and the 
Chern-Simons (CS) part S$ s are expressed as 



S 



DBI 



N C G 



YM 



d A xdz 



K 2 (z)M 2 KK Tv[T^] 



S 



as 



w 5 (A). 



(30) 
(31) 



where the rescaled 't Hooft coupling constant is defined 
as Gym = A/(1087r 3 ) and the field strength of the 5D 
gauge fielcfl A M is Tmn = dMA N -d N AM-i[AM,AN]- 
Here, K\^.(z) are the metric functions of z constrained 
by the gauge/gravity duality. The gravity enters in the z 
dependence of the YM coupling giving rise to the warp- 
ing of the space. In Eq. ([3"T]) . M 4 and R stand for the 
four-dimensional Minkowski space-time and z-coordinatc 
space, respectively, and w 5 (A) is the CS 5-form written 
as 



w 5 (A)=Tv 



AT 2 



-A A T 



1 

10" 



zA l 



(32) 



3 We use the index M = (p, z) with fi = 0, 1, 2, 3. 



Here, T = dA + iAA is the field strength of the 5D gauge 
field A = A^dx^ +A z dz. It should be noted that the DBI 
part is of 0(X 1 ) while the CS term is of 0(A°) with the 
't Hooft coupling constant A. 

We should stress here that as noted in Ref. Q, the 
structure of the the action ([29]) is shared both by the top- 
down Sakai-Sugimoto model and the bottom-up models 
such as in Refs. [HI, d|| as well as the moose models 
in Ref. fill ], with the difference appearing only in the 
warping factors. This allows us to write down a "mas- 
ter formula" which applies to all holographic models and 
moose construction given appropriate warping factors. 

Now to induce the HLS Lagrangian from the action 
(|30|) . we use the mode expansion of the 5D gauge field 
Am(x, z) and integrate out all the modes except the pseu- 
doscalar and the lowest lying vector mesons, which re- 
duces Am(x, z) to A™j° s (x, z). In the AJx, z) = gauge, 
this implies the following substitution @, 3 O 



A.IIWV'IOO, 



(33) 



where {ipn} are eigenfunctions satisfying the following 
eigenvalue equation obtained from the action (|30l) . 



K^(z)d z [K 2 (z)dMz)] = K1> n (z), 



(34) 



with A n being the n-th eigenvalue (Ao = 0). By sub- 
stituting Eq. ([33} into the action in Eq. ([30]), the HLS 
Lagrangian up to 0(p 4 ) can be obtained. The explicit 
expressions for the LECs we need are derived as [!, H[ 



fl = N C G YM M 2 KK / dzK 2 {z) Vo(^) 



aft = N c G YM Ml K \M), 

- 2 =N c G YM (tf), 
T 

vi = -y 2 = -iV c G YM ((i + Vi -^o) 2 ) , 

% = -tf4 = -N C G YM (fi (1 + ^) 2 ) , 
Vh = 2 % = = -2N C G YM (V>iVo} j 

v& = - (y 5 + vi) . 

y 7 = 2N C G YM (V-i (1 + Vi) (1 + A - V'o)) , 

z 4 = 2N c G YM (^(l + Vi -■$>)>> 

z B = -2N c G YM (tf(l + il>i)), 

1 
2 



Ml ( + g^i 



-5^0 + ^1 + ^1 + 5 



4 As emphasized in Ref. Q, the procedure of "integrating out" 
adopted here is different from the "naive truncation" that vio- 
lates the chiral invariance. 



;Mi 



(35) 



where Ai is the smallest (non-zero) eigenvalue of the 
eigenvalue equation given in Eq. (|34[) . and () and (()) 
are defined as 



(A) = 
((A)) = 



dzK 1 (z)A(z), 
dzA(z) 



(36) 



for a function A(z). Equation ([33]) provides the master 
formula for the LECs in the HLS Lagrangian induced 
from general hQCD models. Namely, one just needs to 
plug the warping factor and the eigenfunctions of Eq. (|34[) 
into Eq. ([3"5"|) to obtain the values of the LECs. For 
example, Ki(z) = K^ 1 ' 3 (z) and K 2 (z) = K(z) with 
K(z) = 1 + z 2 correspond to the Sakai-Sugimoto model. 

In addition to the general hQCD models, we also con- 
sider the BPS model studied in Refs. [13, H3 which is 
characterized by the flat space-time. In this case, instead 
of solving the eigenvalue equation, the 5D gau ge fi eld is 
expanded in terms of the Hermite function ip n (27[ r28j , 



Ai( z ) 



(-1) 



n-l 



v/(n- 1)!2™- 1 % /^ 



-z 2 /2_ 



-in—1 



dz 



n-l 



(37) 



where n > 1 and ip 1 corresponds to the wave function 
of the lowest lying vector meson. The wave function of 
the Nambu-Goldstone pseudoscalar boson is expressed in 
terms of the error function erf(z), 



ip (z) = erf(z) = -= 



(38) 



In the following calculation, we use the Hermit function 
and the error function as wave functions of the vector 
mode and pseudoscalar mode, respectively. Then, the 
LECs of the HLS Lagrangian are determined by using 
the above tpo(z) and V , i( z ) m to the master formulas in 
Eq. (2D with K(z) = 1. 



B. The a independence 

In the phenomenological analysis, it is well known 
that the HLS parameter a plays an important role 0- 
[j~ol |. With the leading Lagrangian at 0(p 2 ), the choice of 
a = 2 reproduces the Kawarabayashi-Suzuki-Riazzudin- 
Fayyazudin (KSRF) relation and the p meson dominance 
in the pion electromagnetic form factor. At 0(p 4 ), it 
was shown that the quantum correction enhances the in- 
frared value of a, and, therefore, a good description of 
low-energy phenomenology can be achieved with the bare 
value of a being < 2 [IcJ ■ In the holographic approach, on 
the other hand, the parameter a is attributed to the nor- 
malization of the 5D wave function tp^z), which cannot 
be determined from the homogeneous eigenvalue equa- 
tion (|34[) . As a result, it turns out 3 that the physical 
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quantities are independent of the parameter a as far as 
the leading order in N c is concerned!! 

In order to explicitly see that any physical quantities 
calculated with the HLS Lagrangian induced from hQCD 
models are actually independent of the parameter a, we 
start with Eq. (f33|) . We first note that the vector me- 
son mass and the pion decay constant are related by the 
relation, m 2 = ag 2 f 2 , and m p and / w are fixed by their 
experimental values. Therefore, the HLS parameter a 
and the HLS gauge coupling g are connected through 
ag 2 = m 2 / f 2 . Therefore, the g- independence of physical 
quantities is equivalent to their a-indcpendcncc. 

To see the a-independence explicitly, we define 



so that the new function ipi is normalized as 
N C G YM { dzK x (z) \Mz)} 2 = 1. 



(39) 



(40) 



In terms of the normalized wave function ipi( z ) the 5D 
gauge field of Eq. (j3"3")l is written as 



A^(x,z) = v(x)^ ( Z ) + [a M \\(x) + V^x)} 



(41) 



where d^u = (l/g)a^. In terms of the radial wave func- 
tions of the soliton, F(r), W(r), and G(r), we have 



where 



«|| = -«|| = ~G(r) (r x t) 
9 z 



G(r) = — \G(r) + 1 - cos F(r) 
gr 



(42) 
(43) 

(44) 



From the boundary conditions in Eq. (TTS)) . W(r) and 
G(r) satisfy the following boundary conditions: 



W'(0) = 0, W(oo) = 0, 
G(0) = 0, G(oc) = 0. 



(45) 



The coupled equations of motions for F, G and W can 
be obtained by substituting the expression in Eq. (|4"Tj) 
into the action (|29j) and the minimizing the resultant en- 
ergy. This implies that there are no non-trivial boundary 
conditions to determine the absolute sizes of the vector 
meson contributions W(r) and G(r). The normalization 
of the wave functions for the vector mesons W(r) and 
G(r) are fixed from the boundary condition of the pion 



5 When we include the loop corrections which can be regarded as 
a part of l/N c corrections, the a-dependence becomes relevant 
to the physical quantities through the loop corrections. 



contribution F(r) through the coupled equations of mo- 
tion of F(r), W(r), and G(r), which are independent of 
the gauge coupling constant g since they are expressed in 
terms of the normalized xp v This can be restated as fol- 
lows: The LECs of the HLS Lagrangian are determined 
by substituting the expression in Eq. (|41|) into the action 
(|2T)|) , and the LECs are expressed in terms of normalized 
wave functions ip 1 and ip . Then, all the expressions in 
Appendix [A] can be rewritten in terms of F(r), W(r) and 
G(r) without the gauge coupling constant g. As a result, 
the soliton mass is free from the ambiguity of the normal- 
ization of tp 1 , so that it is independent of the parameter 
g as we shall see in the next section. 

Similar arguments also applies to the moment of iner- 
tia. Using Eq. (|2"Tj) one has 



(46) 



where 



— 0/ ~o 

CI|| = <2|| 



t ■ n^(r) + t ■ r ft ■ r£ 2 (r) 
<p(r) 



with 



2r 



1 



■(nxfy, 



(47) 



Zi(r) = -Mr)-l + coBF(r)] 
9 

Ur) = - g Mr) + l-cosF{r)] 



(48) 



From Eqs. (TTg]) and (|2"2"j) . the boundary conditions for £i j2 
read 



Ci^O) = 0, 6,2(00) = 0. 



(49) 



Again all the equations in Appendix [B] can be expressed 
in terms of £i j2 , V, G, W and F. A similar argument to 
that made following Eq. (|45|) yields that the moment of 
inertia should also be independent of the parameter a. 



C. The CS term and the u> meson 

Another important point which should be addressed 
here with the view point of gauge/gravity duality is that 
the CS term is responsible for the role of the u meson 
in the Skyrmion structure. This can be seen by decom- 
posing the 5D gauge field A into the SU(2) and U(l) 
components as 



A — J 4 SU(2 ) + -Au(i). 
Substituting this in the action of Eq. (|30)) leads to 



(50) 



S 5 DBI = N c 



d A xdz 



1 



- -K^z) [Tr (F^FH + Tr (i^F^) 
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(z) [Tr (F^F^) + Tr (f f F 2 )] | 



(51) 



where Fmn is the field strength of the SU(2) gauge field 
^su(2) and Fmn stands for that of the U(l) gauge field 
^4y(i)- This explicitly shows that, without the CS term, 
when the 5D model is de-constructed from the 4D model, 
only the kinetic and mass terms of the iso-scalar vector 
meson uj are allowed. This conclusion can be explicitly 
verified by using a specific hQCD model such as the SS 
model. 

As can be read from Appcndcx [Al the contribution 
from the kinetic and the mass terms of the uj meson to 
the soliton mass is 



= 47T 



dr 



\r 2 (ag 2 f 2 W 2 



W 



which gives the equation of motion of W as 



W" = ag 2 f 2 W- -W 
r 



(52) 



(53) 



in the absence of the CS term. By making use of the 
partial integration with the boundary conditions given 
in Eq. (|T5|) . Mjf ol can be calculated as 



M 8 " o] = 47T 



0. 



dr 



r 

2 



ag 2 f 2 W- -W -W")W 
r 



(54) 



because of the equation of motion of Eq. (|53")) . Therefore, 
in the absence of the CS term, the uj field decouples from 
the other fields and M^ ol vanishes. This is consistent with 
the earlier studies on the Skyrmions stabilized by vector 
mesons [20j . As can be seen from the equation of motion 
of W(r) given in Appendex [A] the hWZ terms provide 
the source terms of the uj meson field. Therefore, in the 
absence of the hWZ terms, the uj field decouples and does 
not contribute to the soliton formation. 



D. The effective Skyrme parameter e 

Finally we estimate the Skyrme parameter e in the 
original Skyrme model by integrating out the isovector 
p meson from the HLS. The original Skymre model La- 
grangian reads 

£ sk = ^Tr {dJJd»rf) + ^Tr [d^UU\ d„UU^} 2 , 

(55) 

where the chiral field U is U = £ 2 and the first term is the 
nonlinear sigma model Lagrangian that can be written as 
f 2 Tr(aj_^a^_), where a_|_ is defined as d ± without the 
vector field. In the earlier analyses [29[ with the HLS 
Lagrangian up to 0(p 2 ), it is known that the Skyrme 
term can be obtained from the p meson kinetic energy 
term in the limit of infinite p meson mass. In this case, 



the Skyrme parameter e becomes the p meson coupling, 
so that we have e = g ~ 6, which is close to the empirical 
value e = 5.45 that is determined from the A-N mass 
difference. In the HLS Lagrangian up to 0(p 4 ), however, 
we have additional contributions from the pure 0{p A ) 
terms that lead to the Skyrme term. Explicitly, after 
integrating out the p meson, the effective Lagrangian is 
obtained as 



£chPT = fl Tr [a ±/J a£] 

1 £4 _ V\ - V2 

2g 2 2 



+ 



T^{^x^ a ±uY , (56) 



where [ , ] is the commutator and { , } is the anticom- 
mutator. The second terms is the Skyrme term and we 
can read the effective Skyrme parameter e as 



1 



1 



2/1-2/2 



2e 2 2g 2 



(57) 



Since the gauge/gravity duality implies that y\ = —y%, 
the last term of Eq. ([51))) vanishes. Using Eq. (f5"T)) and 
the analytic expressions for the LECs given in Eq. (|35[) . 
the Skyrme parameter is written as 



2e 



1 N C G YM i h ;2 n2\ 



(58) 



With the experimental values of the two inputs / w and 
Trip, we obtain the Skyrme parameter e as 



e ~ 7.31. 



(59) 



in the SS model, while in the flat space-time case, i.e., 
for the BPS soliton model, we obtain 



10.02. 



(60) 



These values are larger than the empirical value of the 
Skyrme parameter e = 5.45 because of the contributions 
from the y 1 , y 2 , and z 4 terms that are of 0{p i ). 

Since the moment of inertia I is proportional to 1/ e 3 
in the Skyrme model, with a larger value of e, we have a 
smaller moment of inertia, which results in a larger mass 
splitting between the A and the nuclcon as is verified 
numerically in the next Section. 



IV. NUMERICAL RESULTS FOR THE 
SKYRMION 

In this section, we present the results of numerical cal- 
culations on the Skyrmion properties in the framework 
of the HLS discussed in the previous section. The HLS 
Lagrangian up to 0(p A ) in Eq. (J7J which is considered in 
the present calculation contains 17 parameters, namely, 
A, a, g, y l (i = 1,...9), z 4 , z 5 , and c lj2j3 . We deter- 
mine all these LECs through hQCD models, which are 
characterized by the warping factor K(z) and the wave 



8 



functions ipo and ipi- Then all the LECs are obtained 
through the master formulas given in Eq. ([55)1 , which 
contain the mass scale Mkk, the 't Hooft coupling A (or 
Gym); and the integrals of the warping factor K(z) and 
the wave functions tp and ip-^. In the present work, we 
consider two hQCD models, the SS model and the BPS 
model. 

In hQCD models, Mkk and Gym are free parameters. 
In the present work, we fix them by using the empirical 
values of /„• and m p : 

m p = 775.49 MeV, 
U = 92.4 MeV. (61) 

Then we have complete information to calculate all LECs 
of the HLS Lagrangian through the master formulas. 
Note, however, that the master formulas can determine 
only the product of ag 2 = m 2 / f 2 , and, therefore, a or 
g remains unfixed. However, as discussed in the previ- 
ous section, the physical quantities are independent of 
a (or g). To be specific, we will first work with a = 2, 
which is widely used in the model of the 0(p 2 ) HLS La- 
grangian [T3| , and then examine how each components 
of the soliton mass and the moment of inertia behave as 
the value of the HLS parameter a is varied. This verifies 
numerically how the a independence comes about. 

In the present work, we consider three versions of the 
HLS model induced from each hQCD model. The first 
version is the model that includes the pion, p meson, and 
u) meson. The second one is the model without the hWZ 
terms, i.e., the model that includes the pion and the p 
meson. The third one is obtained by integrating out the p 
meson in the second version of the model. Therefore, this 
corresponds to the original Skyrmc model but with the 
Skyrmc parameter determined by the hQCD model. In 
this section, we will examine the three versions of the SS 
model and of the BPS model. The obtained results will 
be compared with those of the 0(p 2 ) models, such as the 
p-stabilized model of Ref. [l9| and "the minimal model" 
of Ref. (20j that includes the ui meson in a minimal way. 



A. Skyrmion in the HLS induced from the 
Sakai-Sugimoto model 

In this subsection, we first consider the Sakai-Sugimoto 
model [30l . l3l| to determine the LECs of the HLS La- 
grangian. This model is characterized by the following 
warping factor: 

K 1 ( z ) = (i + z 2 y 1/3 , 

K 2 {z) = l + z 2 . (62) 

Since Mkk and Gym are determined by f„ and m p , all 
LECs except a or g can be determined. This will be 
called HLSi model @. As we discussed above, we take 
the commonly used value a = 2 as a typical example and 
then we will test the results by varying the value of a. 




F(r) 

G(r) 

W(r) 



12 3 

r(fm) 

FIG. 1. The soliton wave functions obtained in the 
BLSi(n,p,uj) model with a = 2. F(r), G(r), and W[r) are 
given by the solid, dashed, and dot-dashed lines, respectively. 
Here, F(r) and G(r) are dimensionless, but W(r) is in unit 
of 1/fm. 



The values of the LECs obtained with a = 2 are given in 
the first row of Table. |U 

Equipped with the numerical values of the LECs given 
in Table U the equations of motion for the soliton wave 
function and for the soliton excitations can be solved nu- 
merically, which allows us to calculate the soliton mass 
and the moment of inertia. The main results of the 
present work are summarized in the columns of HLSi 
of Table HU The results of two models with the HLS of 
0(p 2 ) are also presented for comparison. 

The obtained soliton wave functions for the 
HLSi(7r,/9, w) model are shown in Fig. [TJ and Fig. [5] 
This model results in the soliton mass M so \ w 1184 MeV 
and the moment of inertia X 0.661 fm that leads 
to Am ~ 448 MeV. These numbers should be com- 
pared with the empirical values, M so \ = 867 MeV 
and Am ~ 292 MeV. Compared with the widely used 
"minimal model" of the HLS up to <D(p 2 ) [H [H |H, 
this shows that the HLSi(7r, p, uj) model improves the 
soliton mass. 

We then consider the HLSi(7r, p) model that is con- 
structed from the HLS Lagrangian without the hWZ 
terms. In other words, we set c x = c 2 = c 3 = and 
remove the u meson mass term and its kinetic energy 
term in the HLSi(7r, p, w) model to obtain the HLSi(7r, p) 
model. Therefore, this model is very similar to the model 
studied in Refs. 0, In addition to the soliton 

mass, however, we also calculate the moment of iner- 
tia which was not given in Refs. [Ill H3j. And then we 
finally consider the HLSi(7r) model that is defined with 
the Lagrangian (|56[) with the Skyrme parameter given in 
Eq. ([5"9"]) . All the results are summarized in Table ILTl and 
here are several comments made in order. 

1. As claimed in the literature [H H3, [H, H3], we 
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TABLE I. Low energy constants of the HLS Lagrangian at 0(p 4 ) with a — 2. 



Model 


2/i 


2/3 


2/5 


2/6 


Z 4 


Z5 


Cl 


C2 


C3 


SS model 


-0.001096 


-0.002830 


-0.015917 


+0.013712 


0.010795 


-0.007325 


+0.381653 


-0.129602 


0.767374 


BPS model 


-0.071910 


-0.153511 


-0.012286 


-0.196545 


0.090338 


-0.130778 


-0.206992 


+3.031734 


1.470210 



TABLE II. Skyrmion mass and size calculated in the HLS with the SS and BPS models with a = 2. The soliton mass M so i and 
the A-N mass difference Am are in unit of MeV while y/ (r 2 ) w and yj (r 2 ) E are in unit of fm. The column of 0(p 2 ) + oj (j _B m is 
"the minimal model" of Ref. [2(| and that of 0(p 2 ) corresponds to the model of Ref. [IjJ. See the text for more details. 





HLSl(7T,p,w) 


HLSi (n, p) 


HLSi(tt) 


BPS(tt, P ,uj) 


BPS(7r,p) 


BPS(tt) 


0(p 2 ) + uj^B" [20J 


0(p 2 ) [19] 


M so l 


1184 


834 


922 


1162 


577 


672 


1407 


1026 


Am 


448 


1707 


1014 


456 


4541 


2613 


259 


1131 




0.433 


0.247 


0.309 


0.415 


0.164 


0.225 


0.540 


0.278 




0.608 


0.371 


0.417 


0.598 


0.271 


0.306 


0.725 


0.422 



1.5r 




- Ur) 
(p(r) x 10 



-0.5- ' 



r(fm) 

FIG. 2. The excitations of the soliton profile, £ x (r) (solid line), 
£2(7) (dashed line), and y(r) (dot-dashed line) calculated in 
the HLSi(7r, p, u) model with a — 2. Here, £i(r) and ^2( r ) 
are dimensionless, while ip(r) is in unit of fm. 




HLS^TI) 

HLSjCn.p) 

HLS^p.cfl) 



r(fm) 



FIG. 3. Comparison of the soliton wave functions F(r) 
and G(r) in the three models, HLSi(7r), HLSi(7T, p), and 
HLSi(7r, p, w), which are represented by the solid line, dashed 
lines, and dotted lines, respectively. 



found that the inclusion of the p meson reduces 
the soliton mass. In the present work, the soli- 
ton mass reduces from 922 MeV in the HLSi(7r) 
to 834 MeV in the HLSi(7r, p), which confirms the 
claim that the inclusion of the p meson makes the 
Skyrmion closer to the BPS soliton. However, when 
we include the uj meson, the soliton mass increases 
to 1184 MeV. This is in contrast to the naive ex- 
pectation that including more vector mesons would 
decrease the soliton mass. Since the uj meson inter- 
acts with the other mesons through the hWZ terms, 
this observation shows the importance of the hWZ 
terms in the Skyrmion phenomenology. The role of 
the uj meson in the Skyrmion mass and size can also 
be verified by comparing the soliton wave functions 
shown in Fig. [3J This figure shows that the p me- 
son shrinks the soliton wave functions, which can 
be seen by comparing the results from the HLSi(7r) 



and the HLSi(7r, p) models. However, as can be 
seen by the dotted lines, inclusion of the oj meson 
expands the wave functions. All these behaviors 
can be found in the rms sizes */ (r 2 ) w and */ (r 2 ) E 
in Table |TIJ Therefore, we conclude that the p me- 
son decreases the soliton mass while the u> meson 
increases it. 

2. In the moment of inertia, or in the A-N mass dif- 
ference Am, through the collective quantization, 
the role of the p and u> mesons are the opposite to 
the case of the soliton mass. The mass difference 
Am increases by the inclusion of the p meson, i.e., 
from 1014 MeV in the HLSi(tt) to 1707 MeV in 
the HLSi(7r,p), which worsens the situation phe- 
nomcnologically. Furthermore, in the nuclcon and 
A masses, the rotational energy at 0(1/ N c ) is even 
larger than the soliton mass that is of 0(N C ). This 
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FIG. 4. Same as Fig.[3]but for ^(r) and £ 2 (r). 



raises a serious problem to the validity of the col- 
lective quantization method in these models. How- 
ever, inclusion of the tu meson reduces Am and 
the rotational energy appreciably. The compar- 
ison of the soliton wave functions obtained from 
the 0(1/ N c ) rotational energy can be found in 
Fig.|H This shows that £ x 2 (r) are expanded in the 
HLSi (-7T, p, a;) model than in the HLSi (tt, p) model. 
This leads to a larger value for the moment of in- 
ertia in the HLSi(-7r, p, uj) model, which leads to a 
smaller Am through the relation given in Eq. ([25]) . 

The contributions from each term of the La- 
grangian ([7]) to the soliton mass and the moment 
of inertia are also analyzed as functions of the HLS 
parameter a. The results are summarized in Fig. [5] 
The contributions from the 0(p 2 ) terms, 0(p 4 ) 
terms, and the hWZ terms are represented by the 
dotted, dashed, and dot-dashed lines, respectively, 
while the solid lines are their sums. We first verify 
that the contribution from the 0(p 2 ) terms to the 
soliton mass increases with a. On the contrary, the 
contribution from the 0(p 4 ) terms has a negative 
slope with a and its magnitude is smaller than the 
0(p 2 ) termsH which shows that this order counting 
is reasonable in the Skyrmion mass and size. How- 
ever, the contribution from the hWZ terms that are 
connected to the uj meson is highly nontrivial. In 
particular, its contribution is stable as a becomes 
smaller while the 0(p 2 ) contribution decreases. As 
a result, when a — > 1, which corresponds to the 
value in nuclear medium^ [Toj . the contribution of 
the hWZ terms is close to that of the 0(p 2 ) terms. 



6 The magnitude of the 0(p 4 ) contribution is about 15% of the 
0(p 2 ) contribution at a = 2. 

7 1/N C corrections are expected to be highly important in medium, 
so this feature should be taken with caution. 





' 1 

0(p 2 ) 


' 1 
hWZ 
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-- 0(p 4 ) 
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FIG. 5. Dependence of the soliton mass and the moment of 
inertia on the HLS parameter a in the HLSi(7r, p, oj) model. 
Dotted, dashed, and dot-dashed lines are the contributions 
from the 0(p 2 ), 0(p 4 ), and hWZ terms. Solid lines are their 
sums. 



This evidently shows that the role of the w me- 
son may be even more addressed in nuclear matter. 
Therefore, it is highly desirable to investigate the 
role of the u> meson in more detail in Skyrmion 
matter. Our analysis shows that the three com- 
ponents of the Skyrmion mass represented by the 
dotted, dashed, and dot-dashed lines in Fig. [5] have 
very different behavior with a, but their sum is in- 
dependent of the parameter a. Similar conclusions 
can be drawn from the decomposition of the mo- 
ment of inertia as well. Here, we found a slight 
dependence on a, which might be related to the ap- 
proximate method adopted by the collective quanti- 
zation method. More detailed studies on the quan- 
tization method is, therefore, desirable. 

All these observations show the importance of the uj 
meson in Skyrmions. The uj meson increases the soli- 
ton mass and decreases the moment of inertia, which is 
exactly the opposite to the role of the p meson. Further- 
more, only when the uj meson is included, the rotational 
energy is smaller than the soliton mass and thus the stan- 
dard collective quantization can be justified. 



B. Skyrmion in the HLS induced from the BPS 
model 



It was claimed in Rcfs. [27|,|28( that the BPS Skyrmion, 
i.e., the soliton in the flat space 5D YM action, has the 
potentially important feature in the Skyrmion structure. 
In this subsection, we determine the LECs of the HLS 
with the flat space 5D YM action, which we call the BPS 
model. To investigate the warping factor effect we con- 
sider a gauge theory in flat 5D Minkowski space-time. In 
the sense of large 't Hooft parameter A expansion, the flat 
space-time means that we are going to consider the O(A) 
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FIG. 6. Soliton wave functions F(r) and G(r) in BPS(tt), 
BPS(7r,p), and BPS(7r,p,u;) models, which are represented 
by the solid line, dashed lines, and dotted lines, respectively. 
W(r) is in unit of 1/fm. 
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FIG. 7. Comparison of the soliton wave functions F(r) 
and G(r) in the three models, BPS(7r), BPS(-7r,p), and 
BPS(7r, p, w), which are represented by the solid line, dashed 
lines, and dotted lines, respectively. W(r) is in unit of 1/fm. 



\ \ 

\ :•. ¥>(r) 

.v 

• » ^!fn 

/ 

/ 

r / 
/ .• 
/ / 
/ / 



terms since the warping effect is at O(A ). In flat space- 
time, if the infinite tower of the KK mode is included and 
the CS term is turned off, the Skyrmion solution becomes 
the so-called BPS Skyrmion. In Refs. [HHH the infinite 
tower is truncated to include low-lying isovector vector 
mesons. Here, we include the CS term to investigate the 
lo meson effect. The flat space-time is defined by 

K x {z) = K 2 (z) = 1. (63) 

As in the previous subsection, we consider the three mod- 
els, namely, BPS(tt), BPS(tt,p), and BPS(?r, p, to). The 
obtained LECs and the soliton properties are also pre- 
sented in Table U and [Til respectively. The soliton wave 
functions in the BPS models are shown in Figs. [6] and [7] 
as well. 

Wc find that the role of the vector mesons is similar to 
the case of the HLSi model. Namely, the p meson shrinks 
the soliton while the lo meson expands it. Also, without 
the lo meson, the rotational energy, which is 0(1/N C ) in 
the baryon mass, is much larger than the soliton mass 
that is at 0(N C ). This again raises a serious questions 
on the validity of the collective rotation in the absence 
of the uj meson. 

It is also interesting to note that the soliton mass and 
the moment inertia obtained in the BPS(7r, p, oj) model 
are similar to those of the HLSi(7r,/9, w) model, while 
the difference of the corresponding results in the models 
without the u> meson is quite noticeable. Although the 
HLSi(7T, p, ui) and the BPS(-7r, p, uj) models give similar 
results for the soliton mass and the moment of inertia, 
the obtained soliton wave functions are very different. 
To understand this coincidence, we compare the soliton 
wave functions in these two models in Figs. [8] and which 
clearly show the difference between the two models, es- 
pecially in the wave functions of the p meson. The sign 
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FIG. 8. Comparison of the sol wave functions F(r), G(r), 
and W(r) calculated in the HLSi(7r, p, ui) and BPS(7r,p, lo) 
models. 



difference of the wave functions of is due to the dif- 
ferent sign in the source terms, i.e., the signs of c^s in 
Table [TT] and it can be seen that their magnitudes are 
similar in Figs. |8] and GO 

The difference between the HLSi(7r, p, lo) and the 
BPS(7r, p, lu) models can be easily seen in the breakdown 
of the soliton mass and the moment of inertia. Shown in 
Fig. [TU] are the contribution from 0(p 2 ), 0(p 4 ), and the 
hWZ terms to these physical quantities. We first found 
that the dependence of these quantities on a is very sim- 
ilar to that of the HLSi model. On the contrary to the 
HLSi(7r, p, uj) model, however, the 0(p 4 ) contribution is 
as large as 50% of that of the 0(p 2 ) terms at a = 2. 
Therefore, although the obtained soliton mass and the 
moment of inertia have similar values in both models, 
their breakdown clearly shows their difference. Since the 
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FIG. 10. Same as in Fig. [5] but in the BPS(7r, p, uj) model. 



0(p 4 ) contribution in the BPS model is not suppressed 
enough, it would be interesting to study the contribu- 
tions from the higher order terms to see the convergence 
of the BPS model. 



V. SUMMARY AND DISCUSSION 

In this paper, we have investigated the role of vector 
mesons in the Skyrmion properties based on the HLS 
Lagrangian up to 0(p 4 ). All the LECs of the HLS La- 
grangian could be determined sclf-consistently from a 
class of holographic QCD models by making use of the 
general master formulas p5[) . In the present work, we 
considered two hQCD models, namely, the SS model and 
the BPS model. Equipped with the LECs of the HLS 
Lagrangian determined in this way, we have computed 
the Skyrmion properties and compared the results with 
those of the models of the 0(p 2 ) HLS. 

The results summarized in Table [TT] clearly show the 
important role of vector mesons, in particular, that of u j 



vector meson. As claimed in the literature [16ll27l .l28l. l3^ 



we confirmed that the inclusion of the p meson reduces 
the size of the soliton and decreases the soliton mass so 
that the obtained Skyrmion mass becomes closer to the 
Bogomol'nyi bound. It also decreases the moment of in- 
ertia, which leads to a larger value of the A-N mass split- 
ting. However, in the model of pion only or of pion and 
p meson, the obtained moment of inertia is very small 
so that the rotational energy of 0(1/N C ) becomes even 
larger than the soliton mass of 0(N C ). This strongly 
raises the question on the validity of the collective rota- 
tion in these models[f| 

Then we have studied the role of the uj meson that cou- 
ples to the 7r and p mesons through the hWZ terms in the 
HLS, which are induced from the CS term in hQCD mod- 
els. Contrary to the role of the p meson, the uj meson 
inflates the soliton size and increases the soliton mass, 
while it reduces the rotational energy by increasing the 
moment of inertia. Only when the uj meson is included, 
the rotational energy is smaller than the soliton mass, 
which validates the application of the collective quanti- 
zation. This shows that the uj meson has an important 
role in phenomenology as well. 

We further confirmed that the obtained Skyrmion 
properties are independent of the HLS parameter a. This 
is the consequence of the generic properties of the hQCD 
models related to the normalization of the 5D wave func- 
tions. Given that this relation holds in the large N c limit, 
the a-independence will necessarily break down in dense 
medium where 1/N C corrections are expected to be im- 
portant. 

It should be pointed out that, lessons from nuclear 
physics indicates that, the uj meson is responsible to the 
repulsive interaction which prevents the nuclei from col- 
lapsing and the sigma meson (not the fourth component 
of the chiral four-vector in sigma models but a scalar 
meson relating to the Casimir effects) causes attractive 
interaction and the near cancellation of the these two in- 
teractions gives the small binding energy of nuclear mat- 
ter of ~ 16 MeV. However, it is very difficult to treat this 
Casimir effect given that we have a non-renormalizablc 
theory. The present rough estimation in the Skyrmion 
model obtains the Casimir contribution is ~ 500 MeV 
which goes to the right direction of with a correct order 
of magnitude [36[ . 

Since our study shows the importance of the u> me- 
son in the Skyrmion structure, it is natural to investi- 
gate its effects in the Skyrmion matter. The approach 
adopted in the present paper based on the HLS with 
self-consistently determined LECs from hQCD can be ex- 
tended to the study on dense matter. This can be done 
by constructing the Skyrmions crystal lattice to deter- 
mine the critical density at which a Skyrmion (or an in- 
stanton) transforms into two half- Skyrmions [33, [H| (or 
half-instantons/dyons [3i|). This will be important in 



See, for example, Ref. 1351 for the validity of the collective quan- 
tization in the Skyrme model. 
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understanding of the equation of state for compact-star 
matter as shown in Ref. [4(|. As suggested in Ref. Q, 
a reliable treatment will require low-mass scalar degrees 
of freedom which will figure at subleading order in N c . 
Such work is in progress and will be reported elsewhere. 
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Appendix A: The soliton mass and the equations of motion for F(r), W(r), and G(r) 



Using the wave functions defined in Eqs. (fTBj) and (| 1 T[) and the Lagrangian in Eq. (|7|), the soliton mass in the HLS 
up to 0(p 4 ) is obtained as 

(Al) 



M sol = 4tt / dr [M ( a)(r) + M (4) (r) + M anom (r)] , 
where M( 2 ), Mu\, and M anom are from £(2), £(4)j, + £(4) z , and £ Mom , respectively. Their explicit forms are 
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+ (j/5 + 2/g) 

Ai = z 4 sin F, 
A 3 = 1, 
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The equation of motion of W reads 



W" 



W + ag 2 f 2 W + (y 3 + y A )g 2 W 



g 2 W 2 



G + 2 sin 2 



F 



(2/5 + %) 



g 2 W 



F + — sin 2 F 



(G 2 + 2G + 2 sin 2 F) F' + 2 cos F G + 2 sin 2 - F' + 4 sin FG' 



This evidently shows that the hWZ terms, i.e, the c t terms, are the source terms of W{r). 



Appendix B: Moment of inertia and equations of motion of the excited fields 

When the collective rotation is introduced, the Lagrangian can be written as 

L = -AI sol +XTr(AA^ 
where the moment of inertia X is summarized as 

X = 4tt dr [J( 2 )(r) + J (4 )(r) + JanomM] 



The contributions from £(2), £(4)y + £(4) z , and £ anom , arc represented by I( 2 )(r), F( 4 )(r), and I a nomO 
We further write 
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The hWZ terms give 
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Then the equations of motion for fj, f 2 j and <p are obtained as 
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